Parameter-free density functional for the correlation energy in two dimensions 
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Accurate treatment of the electronic correlation in inhomogeneous electronic systems, combined 
with the ability to capture the correlation energy of the homogeneous electron gas, allows to reach 
high predictive power in the application of density-functional theory. For two-dimensional systems 
we can achieve this goal by generalizing our previous approximation [Phys. Rev. B 79, 085316 
(2009)] to a parameter-free form, which reproduces the correlation energy of the homogeneous 
gas while preserving the ability to deal with inhomogeneous systems. The resulting functional is 
shown to be very accurate for finite systems with an arbitrary number of electrons with respect to 
numerically exact reference data. 

PACS numbers: 73.21. La, 71. 15. Mb 



I. INTRODUCTION 

With the present technology, electron gas can be con- 
fined in various ways to create nanoscale devices of 
lower dimensionality. In particular, the field of two- 
dimensional (2D) physics has grown rapidly alongside 
the development of electronic (quasi-) 2D devices such as 
quantum Hall bars and point contacts, and semiconduc- 
tor quantum dots (QDs). In the modeling of QDs, the 
most common and usually valid approach is to consider 
a purely 2D Hamiltonian with a standard, (not screened 
nor softened) Coulomb interaction and effective values 
for the electron mass and dielectric constant, which ac- 
count for the surrounding semiconducting material such 
as GaAsi 

QDs are studied theoretically in various ways including 
analytic methods, exact diagonalization 2 (ED), variants 
of quantum Monte Carlc"^ (QMC) techniques, Hartree- 
Fock methods 6 , and density-functional theory 7 ^, (DFT). 
The applicability of DFT crucially depends on the ap- 
proximation for the exchange-correlation energy func- 
tional. In spite of recent advances in DFT tailored for 
strongly correlated electrons^ or in the development of 
2D functionals 11 ^— beyond the commonly used 2D local 
spin-density approximatio n 18 ' 19 (LSDA), there is still a 
long path ahead to reach the predictability and efficiency 
that DFT has in quantum chemistry. 

In the present work we focus on the correlation energy 
of inhomogeneous 2D systems within DFT. In particu- 
lar, we consider a generalization of the 2D functional 
developed in Ref. Il4 , which was based on correlation- 
hole modeling similar to that of Beck o 20 ' 21 in 3D systems. 
In this functional several exact constraints are satisfied, 
and the electron spins and currents are incorporated in a 
natural way allowing to deal with spin-polarized and/or 
current-carrying states. However, the functional depends 
on two constants which enter the estimation of the char- 
acteristic size of the correlation hole (see below), and 
they were chosen semi- empirically by fitting the corre- 



lation energy to particular finite systems. Although a 
good performance was obtained, the desired tendency to 
work for an arbitrary number of particles was missed. 
Here, we show how this limitation can be overcome by 
transforming the above-mentioned arbitrary constants to 
non- arbitrary Junctionals of the particle density. This 
is achieved, as explained in detail below, by enforcing 
the functional to reproduce the correlation energies of 
the homogeneous 2D electron gas (2DEG). As a result, 
not only the correlation energies of few-electron QDs are 
reproduced very accurately, broadly outperforming the 
standard LSDA approximation, but more importantly 
the favorable tendency of increasing accuracy with the 
number of electrons is obtained as well. This makes the 
functional a predictive tool to calculate correlation ener- 
gies in realistic 2D systems. 

The rest of the paper is organized as follows. In Sec.lTTl 
we briefly review the framework behind the functional, 
which is similar to that in Ref. 14 In Sec. IIIII we de- 
termine the coefficients for the characteristic sizes of the 
correlation holes by fitting to the correlation energies in 
the 2DEG, both in the fully-polarized and unpolarized 
situations, respectively. In Sec|lV]we apply the resulting 
expressions to QD systems for which accurate reference 
data is available. Conclusions are given in Sec. [V] 



II. MODELING THE CORRELATION HOLE 

We start with the formal expression for the correla- 
tion energy which can be written in Hartree atomic units 
(a.u.) as 



(1) 

where h% a (ri,r2) is the correlation- hole function. In 
Ref. [13 we considered its cylindrical average h° a ' (r = 
rx, s = |r2 — ri|) and constructed a model satisfying the 
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(i) exact normalization of the spin-dependent correlation- 
hole functions (sum rule); (ii) correct short-range behav- 
ior for s — > 0, obtainable from the cusp conditions for the 
2D electronic wave function; and (iii) proper decay in the 
limit s — > oo, for which we used a Gaussian approxima- 
tion. For the correlation-hole potentials 



XJ° a (r) = 2tt / ds h™ (r, s) 



(2) 



we obtained the expressions 
1 fi 

Ur(r) = — (8-37r)^(r)z 2 CT (r) 

817T 

2z aa (r) - 3 In f ~2 CTCT (r) + 1 



(3) 



and 



Ur(r) = (2-^W(r) 

x [2^ ff (r)-ln(2^(r) + l)] (4) 

for the same- and opposite-spin cases, oo' = oo and 
oo' = oo, respectively. The spin-pair components of the 
correlation energy can be calculated from 



(5) 




so that the total correlation energy is given by 

E c [ Pt , Pl ]=Ep + E^ + 2E^ (6) 

with the condition E^ = . 

Let us next examine the ingredients of Eqs. (J3j) and 
©. First, in Eq. © we have 

1 (V/y) 2 _ gA 

4 Pa p a J 

containing two quantities that depend on the oc- 
cupied Kohn-Sham orbitals: t„ = J2k=i I^V'fc.trl 2 
is (double) the kinetic-energy density, and j P)tr = 

^ the spin- 

dependent paramagnetic current density. It can be easily 
seen that D a (r) vanishes for all single-particle (N = 1) 
systems. Hence, £^ CT vanishes as well and the functional 
is correctly self-interaction free for N = 1 (in contrast 
with, e.g., the LSDA). 

As we have mentioned above, important quantities in 
Eqs. © and © are the characteristic sizes of the correla- 
tion holes z aa r(r). They are assumed to be proportional 
to the sizes of the corresponding exchange holes, i.e., 

z aa (r) := 2<w|C£(r)|- 1 , (8) 
z ag (v) := c^O^WI^ + I^WI" 1 ], (9) 

where U£ is the exchange-hole potential 10,22 for spin o. 
The idea behind the assumption is the following: the 
smaller the Fermi hole around each electron is, the more 
tightly the electrons are screened, and therefore they are 
expected to be correlated much less. This coarse picture 
of the real situation has been found to work very well in 
practice^a^S 



III. DETERMINATION OF COEFFICIENTS c atT , 

In Ref. [IH the coefficients c aa and c aa were constants 
determined by fitting the total correlation energy of a 
set of parabolic QDs. In the present work, these spin- 
dependent coefficients are expressed as functionals of 
the particle density. We achieve this goal in two steps. 
First, we determine the values of c aa i that yield the ex- 
act correlation energy density^ of the polarized (un- 
polarized) 2DEG, denoted below as e 2DEG [r s ,C = 1] 
(e 2DEG [r s , ( — 0]), where ( is the spin polarization. This 
allows us to write 

Coo' -> Coo'[r s }- (10) 

where r s = 1/ ^/ifp and p = p^ + pi is the total density of 
the 2DEG. Then as a second step, when the functional 
is applied to an inhomogeneous system, we express r s in 
terms of the local particle density, i.e., for each point in 
space. This is nothing else but a local-density approxi- 
mation for the coefficients c aa > . Otherwise, the correla- 
tion functional has the form given in Eqs. ([2][5]). Hence, it 
still satisfies all the exact constraints listed above and the 
overall expression is of the form of a current-dependent 
meta-generalized-gradient approximation (meta-GGA). 



A. Fully-polarized case 

Now we demonstrate in detail the procedure determin- 
ing c a<7 i — > c aa f[r s ]. First we focus on the fully-polarized 
case with oo' = oo, where a =j" or o —\,, and thus £ = 1. 
In this case p = pf or p = p±. Comparing Eq. ([5]) with 
the LSDA expression, 



E, 



LSDA 



= dvp{v)el^[r s ,( = l} 



(11) 



directly yields = 2e 2DEG [r s ,C = 1]. Similarly 

for the exchange-hole potential in Eq. © (see Ref. [lj 
for its definition) we find U a x = 2e 2DEG [r s ,C = 1] = 
-16/(37rr s ). From Eq. © we get z aa = 2c aa \U^\- 1 = 
3<KC aa r s /8, and Eq. © for the 2DEG gives D a = vr^Vj 4 . 
Collecting these results to Eq. © leads to 



2DEG 



[r.,C = l] = 
-3tt) c\ a 
24 r? 



7T 

1 Caa ' 



In 



(12) 



This expression can be solved numerically for c ff<r [r s ] by 
using the parametrized QMC result for the correlation 
energy of the 2DEG?iS The result is shown in the upper 
panel of Fig. [T]for realistic densities (0 < r s < 20). The 
data could be easily tabulated or parametrized for con- 
venient use of the functional. For this range of densities 
we may use approximate parametrization of the form 



a log(r s ) + I3r] 



(13) 
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FIG. 1: Upper panel: Optimal coefficients c CT(T for locally re- 
producing the correlation energy densities of the homogeneous 
two-dimensional electron gas. The inset show the detailed be- 
havior in the high-density limit r s — ► 0. Lower panel: Same 
for Caa- Solid and dashed lines correspond to two difference 
approximations for the polarized component of the correla- 
tion energy in the unpolarized electron gas (see text). The 
inset corresponds to the high-density result in the first ap- 
proximation (APP1). 



with a = -0.1415 1,0 = 1.226 1, and 7 = 0.14499. We 
remind that when using the correlation functional the 
coefficient c crcr [r s ](r) is calculated locally with the local 
density r s (r). 

Let us examine in detail the high-density limit, r s — > 0, 
for the same-spin coefficient. Expansion of the logarithm 
in Eq. {12} leads to 



-►01 = 



768 



7T 2 (8-37T) 

1.2087, 



£ 2DEG r 



o,c = i; 



1/4 



(14) 



where we have used the known limit for the 2DEG, 
ef EG [r s -> 0,C = 1] ~ -0.0390750 (Ref.QJ). This re- 
sult coincides with the numerical one in Fig. Q] (see the up- 
per inset). We point out that the "u-shape" at r s < 0.05 
is not due to numerical inaccuracies in the present work 



or in the QMC parametrization, but it arises from the 
exact properties of the 2DEG correlation energy in the 
high-density limit.— 

In the low-density limit, r s — > 00, we obtain from 
Eq. (fT2]) the leading contribution as 



96 r. 



tt(8 
1.68, 



3tt) 



-e? DEG [r, 



oo,C = l] 



1/3 



(15) 



where we have used the known result for the 2DEG, 
£ 2DEG[ rs _>. = l] _> -0.222/r s (Ref.Il|). From 

the upper panel of Fig. [1] it is evident that this extreme 
low-density limit has been not reached yet for r s ~ 20. 
However, our numerical values at larger r s agree very well 
with the exact limit. 



B. Unpolarized case 

The coefficient c CTff [^s] can be determined in a similar 
fashion by employing the results for the correlation en- 
ergy density of the unpolarized 2DEG, e 2DEG [r s ,C = 0]. 
It should be noted that also the same-spin component 
U° a depending on the spin density p a is needed in the 
calculation. Now, the fitting to the 2DEG raises the 
following complication: the parallel-spin contribution of 
e™ EG [r s , C = 0], which we here mark as e;?° EG [r S) C = 0], 
is nontrivial; in particular, it is important to note that 
this contribution is not the same as e 2DEG [r s , £ = 1] upon 
density scaling. In principle, a fully spin-resolved expres- 
sion is accessible by employing the QMC results for the 
2DEG correlation potential energy and using the virial 
theorem^ For the sake of simplicity, however, we apply 
here two different approximations for the parallel-spin 
component. In the first approximation (APP1) we set 
e™ EG [r s ,C = 0] w el UEG [V2r s ,C = 1], which is simi- 
lar to the form of Stoll et alM- in 3D electron gas. In 
fact, this approximation corresponds to the density scal- 
ing mentioned above. The second possible approxima- 
tion (APP2) is given by e™G[r s , C = 0] « ef EG [r s , £ = 
l]/y/2, which is similar to the 3D version of Perdew and 
Wang.— Even though the first approximation is exact in 
the limit r s — ¥ 0, the latter type of approximation has 
been shown to be more accurate- 7 in 3D, and we may 
thus expect similar tendency in 2D as well. It should be 
noted, however, that deviations from numerically exact 
results have been reported in 3D for both approxima- 
tions, especially at large r s (Ref. |H). 

Going back to the determination of Co-g-^], we consider 
spin-densities p-[ = Pi = p/2 and set a condition U° a — 



2el D a „[r s ,( = 0] with the approximations APP1 and 
APP2 for the r.h.s. as specified above. From Eqs. 
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(|6j) the formula to be solved for cw^s] is given by 
^ 2DEG ^,C = 0] - ef EG [r s ,C = 0] = 



3iry/2 



-In 1 



3tt\/2 



(16) 



The result is shown in the lower part of Fig. [T] for both 
APP1 (solid line) and APP2 (dashed line). Both curves 
for this density range (0 < r s < 20) can be approximated 
with a satisfactory accuracy by a simple parametrized 
formula 



f [r s ] « Sri 



(17) 

S = 



with S = 0.653 58 and f = 0.116 91 for APP1 and 
0.663 25 and f = 0.123 96 for APP2. 

Again, let us consider the high- and low-density limits 
for the obtained expressions. In the high-density limit 
we find a semi-analytic expression 



0] = 



«"[*■. ->o,c = o] 



3^-2) 

- e f EG [rw0,C = 0] 



1/2 



(18) 



Usingthe known limit ef EG [r s -> 0, ( = 0] « -0.192 500 
(Ref.llJ) leads to c^ pi [r s -> 0] 0.55168, which cor- 
responds to the exact result in this limit as mentioned 
above, and c^? P2 [r s -> 0] rj 0.57175. Both values coin- 
cide with the numerical results in Fig. [T] (see the lower 
panel for a detailed view on APP1). 
In the low-density limit Eq. (|16|) yields 



16 r. 



e^[rwoo,C = 0] 



3\Z2(tt - 2) 

-ef EG [r fl -- x., = l.;[ (!<)) 



r s ->• 00 = 



From the knowledge that e;? DEG [r s — > 00, ( = 0] — > 
-0.470/r s (Ref. M), we obtain CiyS [r s -> 00] ~ 1.03. 
Note that the limit value is the same for APP1 and 
APP2, which is obvious due to the low-density decay 
of the correlation energy as oc 1 jr s in the fully-polarized 
2DEG (see the previous section) . The congruence in the 
limit value between APP1 and APP2 is not evident in 
the density range (r s < 20) shown in Fig. [TJ In fact, the 
decay towards the same limit can be seen only at very 
high r s as we have confirmed numerically. 

It is reassuring to note that, to the leading order, both 
high- and low-density limits of our model correlation en- 
ergies obtained through Eqs. (|T2|) and flTBl are in agree- 
ment with the corresponding exact expansions regarding 
their functional dependence with respect to the density 
parameter r s . More specifically, e 2DEG [r s — > 0, £ = 0] 



£ 2DEG r 



0, £ = 1] attain a finite (r s indepen- 



and 

dent) value in the high-density limit, while £c UEG [ 
00, C = 0] and e 2DEG [r s — > 00, £ = 1] decay as r~ 
low-density limit. 



r s ->■ 
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FIG. 2: (color online). Relative errors produced by the 
present functional (circles) and the local spin-density ap- 
proximation (squares) in correlation energies of spin-polarized 
quantum dots with respect to numerically accurate results in 
Ref. |2 (see text). 



IV. TESTING ON INHOMOGENEOUS 
SYSTEMS 

Next we test the functional, i.e., Eqs. (J3]l9j) with the 
coefficients c aa and c aS taken from Eqs. (|13l) and (fT7|) . re- 
spectively, for a set of 2D parabolic quantum dots, where 
the external confining potential is given by v(r) — w 2 r 2 /2 
with a; being the confinement strength. The correla- 
tion energies obtained from Eq. ([5]) arc compared to 
the reference results E r c Qt = E\l\ — Ef^ x , where E\°\ is 
the analytic,^ ED^ or QMC,i total energy, and £ t E ^ x 
is the exact-exchange (EXX) total energy obtained in 
the Krieger-Li-Iafrate approximation 30 with the octopus 
code.' 51 The EXX result is used as input for our func- 
tional, including ~U° as the exact exchange-hole poten- 
tial. 



A. Fully-polarized quantum dots 

Figure [5] shows the relative errors in the correlation en- 
ergies for fully-polarized QDs of N = 2 ... 6 electrons and 
confinement strengths w = 1/4 (solid lines) and u = 1/16 
(dashed lines) . The errors produced by the present func- 
tional and by the LSDA with respect to the ED results 
in Ref. |2j are presented by circles and squares, respec- 
tively. The present functional clearly outperforms the 
LSDA by reducing the error in the correlation energy by 
one order of magnitude on the average (mean absolute 
errors 120% and 11%, respectively). The performance 
of the present functional is stable regardless of N and w, 
whereas the LSDA gains accuracy very slowly as a func- 
tion of N. Since our functional coincides with the LSDA 
in the 2DEG limit by construction, it can be expected 
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TABLE I: Comparison of the correlation energies for unpolar- 
ized parabolic quantum dots. The reference results have been 
obtained from numerically accurate data for the total energy 
(see text). APP1 and APP2 correspond to our functional with 
two different approximations for the parallel-spin component, 
respectively (see text). The last row shows the mean absolute 



error 


(in percentage), A 


= (\Ec - 


E? { \/\E? l \). 




N 




zpref 
■c-tot 




£ C APP1 


jjAPP2 


77LSDA 


2 


1 


3* 


-0.162 


-0.144 


-0.147 


-0.199 


2 


1/4 


0.9324+ 


-0.114 


-0.098 


-0.101 


-0.139 


2 


1/6 


2/3* 


-0.102 


-0.085 


-0.087 


-0.122 


2 


1/16 


0.3031+ 


-0.070 


-0.057 


-0.059 


-0.085 


2 


1/36 


0.1607+ 


-0.049 


-0.040 


-0.041 


-0.061 


6 


1/1.89 2 


7.6001* 


-0.421 


-0.396 


-0.405 


-0.473 


6 


1/4 


6.995+ 


-0.396 


-0.381 


-0.390 


-0.457 


6 


1/16 


2.528+ 


-0.250 


-0.221 


-0.228 


-0.279 


12 


1/1.89 2 


25.636* 


-0.917 


-0.895 


-0.915 


-1.000 


A 








11% 


9% 


18% 



* Analytic solution by Taut from Ref. |29l + CI data from 
Ref. * Diffusion QMC data from Refill 



that the error of the present functional diminishes fur- 
ther at larger TV. This desired tendency is clearly missed 
if a fixed value for c aa is used^ 



B. Unpolarized quantum dots 

In Table Q] we compare the correlation energies for a 
set of unpolarized QDs. Also in this case the present 
functional, with both approximations APP1 and APP2 
for the parallel-spin component, is more accurate than 
the LSDA. As expected, LSDA becomes again more ac- 
curate with N, but this tendency can be found also in 
the present functional, where it is in fact significantly 
stronger than in the LSDA. For example, for N = 2, 
N = 6, and N — 12 with roughly the same w, the LSDA 
yields relative errors of —22%, —12%, and —9% in the 
correlation energies, respectively. When using our func- 
tional with APP1 the corresponding errors are 14 %, 6 %, 
and 2%, respectively, and within APP2 the errors re- 
duce further to 11 %, 1.5 %, and 0.2 %. As can be seen in 
Ref. 0, the functional having a fixed value for c aS is far 
from this accuracy. The better performance of APP2 in 
comparison with APP1 is in line with the results for the 
3D electron gas, where the approximation of Perdew and 



Wang 26 (similar to APP2) is more accurate than the one 
by Stoll et al& (similar to APP1). 

Before concluding, it is natural to ask if the exchange- 
hole potential U£(r) may be treated approximately to 
reduce the burden of the EXX calculation without losing 
the accuracy. Possible choices may be the LSDA expres- 
sion, or the GGA and meta-GGA as given in Refs. ll^ML3l 
Also the goal is to carry out fully self-consistent calcula- 
tions and preferably for a larger variety of 2D systems, 
e.g., for quantum- ring structures. Here we have focused 
solely on parabolic QDs representing, through compari- 
son with experimental dataj^Sr— a valid approximation 
for both vertical and lateral semiconductor QD devices. 
The above outlined tasks clearly deserve extended future 
investigations. 



V. CONCLUSION 

We have developed a spin-dependent parameter-free 
density functional to calculate the correlation energies 
in two-dimensional electron systems. The functional has 
been constructed through physically reasonable modeling 
of the angular-averaged, spin-resolved correlation hole 
functions. The key extension to previous works is en- 
forcing the functional to reproduce the correlation en- 
ergies of the homogeneous two-dimensional electron gas. 
We have shown that this is possible by transforming the 
coefficients - involved in the estimation of the charac- 
teristic sizes of the correlation holes - to functionals of 
the electron density. As a result, we are able to find very 
accurate correlation energies for quantum dots with vary- 
ing confinement strength and number of electrons. Most 
importantly, the significant error reduction as a func- 
tion of the number of electrons makes the present func- 
tional a predictive method to obtain correlation energies 
of systems which are beyond the capabilities of exact- 
diagonalization and quantum Monte Carlo techniques. 
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